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We consider a model describing the one-dimensional confinement of an ex- 
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00 ' citon in a symmetrical, rectangular quantum-well structure and derive upper 

and lower bounds for the binding energy of the exciton. Based on these 

d ' 

S bounds, we study the dependence of E^ on the width of the confining po- 

tential with a higher accuracy than previous reports. For an infinitely deep 

o" 

potential the binding energy varies as expected from 1 Ry at large widths to 
4:Ry at small widths. For a finite potential, but without consideration of a 
mass mismatch or a dielectric mismatch, we substantiate earlier results that 
the binding energy approaches the value 1 Ry for both small and large widths, 
having a characteristic peak for some intermediate size of the slab. Taking the 
mismatch into account, this result will in general no longer be true. For the 
specific case of a Ga\- x Al x As / GaAs / Ga\- x Al x As quantum- well structure, 
however, and in contrast to previous findings, the peak structure is shown to 
survive. 
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I. INTRODUCTION AND STATEMENT OF THE PROBLEM 



The study of electronic and excitonic properties in quantum well structures has been a 
subject of great interest since the pioneering work of Dingle, Wiegmann and HenryS. In view 
of the enormous amount of literature in this field, any list of references must be incomplete. 
We quote the publicationsEHil, which are in part related to this work, and recommend the 
references therein. One of the most appealing features of these systems is the enhancement 
of excitonic effects, for instance the increase of the binding energy and the oscillator strength, 
which may allow the observation of an exciton even up to room temperature. 

The main reason for the large binding energies and oscillator strengths is understood to 
be the quantum confinement of electron and hole in the growth direction of the heterostruc- 
ture; in comparison to a three-dimensional monostructure, the electron-hole correlation is 
increased. Moreover, the mass mismatch as well as the dielectric mismatch in a quantum 
well may even enhance this effect, as was primarily pointed out by Keldyshlll. 

In contrast to the simple qualitative reasonings, a quantitative theoretical description 
of the excitonic enhancement is quite complicated. The reasons are obvious; as translation 
symmetry is broken in the growth direction of the heterostructure, the familiar separa- 
tion of center-of-mass and relative part of the exciton motion is no longer possible. If 
a dielectric mismatch is to be included, the electron-hole potential is no longer a simple 
Coulomb potential. If a mass mismatch exists, the kinetic-energy part of the Hamiltonian 
is no longer isotropic. Turning to real substances such as Ga\- x Al x As / GaAs / Ga\- x Al x As, 
band- structure complications (e.g. valence-band degeneracy) do occur. Moreover, the grow- 
ing process may induce an interface roughness etc.. As a consequence, the spectra of excitons 
in a quantum well are far from being understood on a quantitative scale and different aspects 
have been controversially discussed in the literature. 

In this paper we are concerned with one of the problems which — to the best of our 
knowledge — has not been solved; i.e. the position of the peak of the excitonic binding energy 
as a function of the width of the confining potential. Some introductory comments may be 
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appropriate. For an infinitely deep well, the binding energy is known to vary monotonically 
from 1 Ry to 4 Ry if the width L varies from infinity to zero. Contrasting this case with that 
of a finite well in an otherwise homogeneous material, the binding energy behaves similarly 
for sufficiently large values of L, but qualitatively different for small L. In the latter case, 
the wave function spills increasingly over the interfaces and into the barriers, occupying 
a greater three-dimensional volume as tunneling becomes more and more important. In 
the ultimate limit of zero width, the binding energy will again be 1 Ry. For a finite well 
width the confinement will increase the binding energy Ej, above 1 Ry. Thus we are lead 
to conclude that Eb should have a maximum for some intermediate well width. In fact, 



Sec. IV, we complement their results by lower bounds for the binding energy, which have a 
relative deviation of (at most) 0.2 from the corresponding upper bounds and, furthermore, 
a similar shape. 

A finite well in an otherwise homogeneous material should probably be viewed as a rather 
poor model for a real quantum well. We anticipate, however, that such a conclusion might 
be somewhat pessimistic. The model is well applicable as starting point for a quantitative 
description, but we include mismatches for the masses and the dielectric constants. What 
about the binding energy under these circumstances? The situation is clear in the limiting 
cases; for infinite well width, we start with 1 Ry of the well material, whereas for zero 
well width, we arrive at 1 Ry' of the barrier material. Whether or not a peak appears 
will be sensitive to the value of Ry/Ry'. Andreani and Pasquarelloi performed a study for 
Gai- x Al x As / GaAs / Gai- x Al x As, including the effects of, firstly, a dielectric mismatch and, 
secondly, the valence band degeneracy. They did not find a peak for x > 0.25 and a width 
L > 30 A, even when the band degeneracy is (theoretically) switched off. 

The intention of this paper is to critically reexamine this conclusion. Clearly, our model 
must also resort to simplifications, as was indicated above. For each part of the heterostruc- 
ture, we assume nondegenerate, isotropic and parabolic bands, but include a mass mismatch 
and a dielectric mismatch at the interfaces. The confinement of electron and hole is mim- 
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icked by finite rectangular wells. Our strategy is to produce upper and lower bounds for 
the correct binding energy of the model. Consequently, we can estimate the error of our re- 
sults quantitatively. Inserting the material parameters of Ga\- x Al x As / 'GaAs / ' Ga\_ x Al x As, 
x varying from 0.15 to 0.40, we do find a peak structure. We can thereby disprove the above 
assertions. 

II. FORMULATION OF THE MODEL 

To begin with, we fix the geometry of the quantum well as shown in Fig. [l|. The figure 
shows a cut perpendicular to the y-axis, the well extending from z = —L/2 to z = L/2. The 
remaining space \z\ > L/2 is occupied by the barrier material. The material parameters of 
well and barrier are denoted by unprimed and primed characters, respectively. To achieve a 
compact nomenclature, we define ^-dependent expressions, e.g. 

m(z) := m6(L/2 - \z\) + m'9(\z\ - L/2), (1) 

for the mass of a particle; 8(z) is the familiar step function. As indicated in the figure, 
we assume translational invariance within the x-y plane; consequently, the corresponding 
components of the total momentum are conserved. 

To set up a model Hamiltonian, we proceed as described before. Introducing center-of- 
mass and relative coordinates for the motion in the x-y plane, the corresponding center of 
mass part can be eliminated by projection on the subspace of total momentum zero. We are 
left with the Hamiltonian 

H := (2p, ± (z e ,z h ))~ 1 p 2 1 _ + p z ,e{'2m e (z e ))~ 1 p z , e + p Zt h(2m h (z h ))~ 1 p Zth 

+V e {z e ) + V h {z h ) + V eh (f ± , z e , z h ). (2) 

Here, the indices e, h characterize electron and hole, respectively; z e , z^, p z , e , and p z> h are 
the z-components of position and momentum, m e (z) and rrih(z) the corresponding masses. 
For the relative coordinate in the x-y plane we use the notation f± := r± te — r± th , p± denotes 
the corresponding momentum, and 
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(n±(z e , z h )) 1 := (m± >e (z e )) 1 + (mx,h(zh)) 1 (3) 

is a generalized reduced mass. The confining potentials for electron and hole are V e (z e ) and 
Vh(zh)- These potentials are assumed to be finite rectangular wells of the width L: 

\ 0, if \zi\ < L/2 , 
Vi(zi):=\ (4) 
[ V t , if \zi\>L/2 , 

where i — e,h. Finally, we have to specify V e h(r±, z e , Zh), that is the potential energy of 
electron and hole. To do so, we apply Poisson's equation to find the electrostatic potential of 
one of the point charges (e.g. that of the hole) under the geometrical conditions of Fig. |l|. At 
the interfaces, the familiar continuity conditions of electrodynamics have to be fulfilled. If we 
assume this potential to be known, its product with the electron charge yields the potential 
energy of interest. We note that the required solution of Poisson's equation can conveniently 
be set up as follows. At first, one performs a Fourier transform of the potential with respect 
to the coordinates x and y. One will realize that the remaining differential equation with 
respect to z can be treated directly (the solutions are exponentials). In fact, this is the 
method, which was used by Fomin and Pokatilov0 under more general circumstances. For 
the present geometry, their formulae can be somewhat simplified, leading to an explicit 
result for the potential energy which can be found in a paper of Kumagai and Tagakahara0. 
These authors use the method of image-charges and arrive at a power-series expansion of 
V e h, the expansion parameter being 



£-£' 



r.—,. (5) 



The zero'th order term V^j of the power series is the familiar Coulomb expression 



e 2 

V$ (r ± , z e , z h ) := , (6) 

er 



the dielectric constant being that of the well material and e denotes the electron charge, 



divided by y/^rreo (SI units). The higher-order terms constitute the so-called image potential. 
For further use, we note the g-linear contribution explicitly: 



(r ± , z e , 2 ft ) := + + . (7) 



e 



In this equation, we denned the distances r, r± and the 9-factors as follows: 



rf + (z e - z h ) 



r± ■= \M + (*> + ^ ± L ) 2 > ( 8 ) 

and 

e>:=^ e -i/2) + 0(z fc -L/2), 

6« :=9(-z e -L/2) + 6(-z h -L/2). (9) 

We add three remarks concerning these formulae. 

The reader will notice that all terms of Vj^ are of Coulomb type, the corresponding 
denominators being partially displaced. It is important to realize that this property remains 
valid for all higher-order contributions (again, see Ref. [11]). Consequently, we induce an 
error of order q 2 , if we replace the full expression for V eh by v}f^ + . As q is typically 
of the order 0.1, this simplification is reasonable and will be used in the remainder of this 
paper. 

The second remark is concerned with an ambiguity in the power-series expansion of V e h- 
Because of Eq. (H), we may replace e by e'{l + q)/(l — q). Inserting this equality into Eqs. 
(|) and (|7|), we find terms of zero'th and first order as follows: 

V e f(f ± ,z e ,z h ):=-^-, (10) 

These expansions are clearly equivalent, and will be used to our advantage in the next 
section. 

Finally, we mention that some authors (see, e.g., Refs. P| JTT|JT^ ) include self-energy terms 
in H, which are also due to a dielectric mismatch. Formally these appear if one considers 
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an electron- and a hole-inhomogeneity in Poisson's equation and inserts the total electro- 
static energy into H . The interaction part of this energy is the one we used above (i.e. 
V e h(rx, z e , Zh)), the diagonal parts produce the self-energy terms £(z e ) for the electron (i.e. 
— V e h(0, z e , z e )/2, evaluated for r i = and Zh = z e , the bulk Coulomb singularity being 
subtracted) and £(2^) for the hole. Inspection of Eq. (0) shows that £(z) consists of one- 
dimensional Coulomb potentials, which exhibit a singularity on the interfaces z = +L/2 and 
z = —L/2 and vanish if q vanishes. Summarizing, the confinement potentials are claimed to 
be changed. We doubt that modifications of the confinement potentials can be motivated 
this way and refer to the early quantum-mechanical debate on the hydrogen problem with 
(or without) self-energy corrections (e.g. Tomonag 80). Consequently, we did not include 
such corrections here. We agree, however, with the statement made in Ref. |l| that a cor- 
rect microscopic description of an exciton in a quantum well should contain (well behaved) 
polarization-induced modifications of the rectangular confinement potentials. It should be 
mentioned that if we modify the confinement potentials in the way indicated above, then 
the corrections to the final binding energies presented in Sec. IV only appear in the vicinity 
of the peak being not larger than 0.5 meV. 

It will prove useful to introduce dimensionless variables. For the excitonic system under 
consideration, appropriate units of length and energy are Bohr radius and Rydberg energy. 
Considering the well material, these are 



t 2 -4 

aB '-=^ Ry:= 2l^- (12) 



Returning to Eqs. (0), (|7|), we replace r by a#r, introduce the dimensionless Hamiltonian 

h := H/Ry (13) 

and find the expression 

h = (i±. ^ 2 _ d n ± d d n ± d 



Vx(z e , Zh) dz e m e (z e ) dz e dz h m h (z h ) dz h 

_2 _ 2q f ®> + ®« + i-e> + i-e« \ 



+U e (z e ) + U h (z h ). (14) 

All variables are now dimensionless. The confinement potentials read as follows; 

0, if \zi\ < 1/2, 
[ U u if N > 1/2, 

where we introduced 

l:=~. (16) 

Ry a B 

The theta factors are defined as in Eq. (^) with the exception that L is replaced by I. 

We shall need a second dimensionless version of Hamiltonian (Q) , which is based on the 
Bohr radius a' B and the Rydberg energy Ry' of the barrier material as well as the formulae 
flnp, (O) for the electron- hole potential energy. Introducing 

ti := H/Ry' (17) 



one arrives at 



h> = ^ y 2 d fi' ± d d fi' ± d 



H±(z e , Zh) dz e m e (z e ) dz e dz h m h (z h ) dz h 

_2 _ 2q + + i - + i-e^ \ 

r \ r r_ r + J 

+U' e (z e ) + U' h (z h ). (18) 

The theta factors are defined as in Eq. (Q) with the exception that L has to be replaced 
by V := L/a' b . In analogy, the confinement potentials U' can be taken from Eqs. flTSP and 
(^), if one replaces the unprimed material parameters by primed ones. 



III. UPPER AND LOWER BOUNDS FOR THE EXCITON BINDING ENERGY 

Because of the absence of translation symmetry, the Hamiltonian H (or h, h' , respec- 
tively) cannot be treated analytically. Nevertheless, we can fix quantitative properties of 
the binding energy Ef, by means of rigorous bounds, which will be derived in the following 
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sections. Our strategy will be to discuss firstly the ground-state energy Eq. The reason is 
that we can directly prove bounds for E Q . The binding energy Eb, in turn, is related to E 
by the equation 

Eb = E cont — E , (19) 

where E cont denotes the energy of the continuum edge. In our case, E cont is the sum of the 
lowest well energies of electron and hole, any correlation being neglected. For further use 
we note the explicit formulae 

E con t/ Ry = e , e + e ,h, (20) 
where e 0ji {i = e, h) is the ground-state eigenvalue of the one- dimensional well Hamiltonian 

^-=-# J ^# + ^W' (21) 
oz mAz) oz 



part of the total Hamiltonian (fn|). The eigenvalue e 0) j is implicitly defined by the equation 



eo,i = ^arcsm =. (22) 

m ^ yjl + {m'Jrm - l)e , t /U t 

Due to Eq. (0), an upper (lower) bound for E yields a lower (upper) bound for E b . 



A. Lower bounds for the binding energy 

We use a variational approach, to derive an upper bound on E . The trial wave function 

is 

*«,A(rL, z e , z h ) := ® e (z e )$ h (z h )e- a Vri+^\ (23) 

where z := z e — Zh and $j (i = e, h) is the ground-state eigenfunction of hi (see Eq. (pTj) ) . 
Clearly, the factors $ e and ^ serve to incorporate the single-particle well behavior. 

Choosing an exponential envelope in Eq. (^), we account for two important aspects. 
Firstly, the confinement-induced change of the effective Bohr radius (parameter a), and, 
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secondly, the quenching of the wave function in ^-direction (parameter A) . Both aspects are 
particularly important if one approaches the quasi two-dimensional case, where the width 
of the layer becomes considerably smaller than the Bohr radius of the exciton. We remark 
that the envelope can correctly reproduce the ultimate limits of a free exciton ground-state 
in two or three dimensions. Subsequently, the trial function is quite flexible and we expect 
the variational inequality 

— < mm — L - J ■ (24) 



Ry \a 

to be effective. For the binding energy Ej, we derive 

— > e e + en h — mm rz . (25) 

Note that the minimum of the right-hand side can be calculated for q = 0, as we have 
linearized h with respect to q. 

B. Upper bounds for the binding energy 

In this subsection, we provide a class of lower bounds for Eq. To do so, we assume that 
the barrier masses m' e , m' h are not smaller than the corresponding well masses m e , m^. This 
is the case for Gai- x Al x As/GaAs/Gai- x Al x As, which will be used as an example. Recalling 
formula (|T8| ) for the Hamiltonian under consideration, one verifies by direct inspection that 
the following inequality is true in the sense of operators: 

ti > h', (26) 



where 



m' e dz 2 e m 'h dz l 
-- -2q(- + -)+ U'M + U' h (z h ). (27) 

We stress that the above assumption vn! i > rrii is not crucial. If the contrary was true, we 
would find an analogous inequality with respect to Hamiltonian ([141) . 
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To provide a lower bound to the spectrum of h', we split the Hamiltonian fl27|) into four 
tractable parts. Assuming the corresponding ground-state energies to be known, their sum 
will be a lower bound to Eo/Ry'. Let us consider the following decomposition: 

h' = K + ti h + ti c + h' im , 

h! i = :-(l-Xi)^j^-z + Uftz i ), i = e,h, 
m - ozf 



K = '■ -Z/V ± - x e y e — — - x h y h 



/4 d 2 n' ± d 2 2 

m'dz 2 h ^ h m' h dzl r 



h u ^h 

2 



m' e dz 2 

-^ l -^ b M- 2q (^ + ^)' (28) 

Here, Xi,yi,y are parameters with values in the interval [0, 1] but otherwise arbitrary, which 
can then be chosen to lift the lower bound as much as possible. The exciton ground-state 
energy E fulfills the inequality 

E Q /Ry'>e' e + e' h + e' c + e' im , (29) 

where the four terms on the right-hand side are the ground-state eigenvalues of the four 
parts of h', defined in Eq. (^). The eigenvalues e- can be derived from Eq. (JZ2"D; removing 
the mass mismatch and rescaling the particle mass appropriately, we find 

e[ = 4(1 - Xl )-^ arcsin 2 sfl - e!jV[. (30) 

Concerning e' c , we transform the corresponding Hamiltonian h' c ; introducing center-of-mass 
and relative coordinates instead of z e and Zh, we can separate the center-of-mass part. 
Performing an appropriate scaling transformation with respect to r± and z = z e — Zh, we 
arrive at the following equations; 

e' = -inf spec J -V 2 - . 2 i =: -e(A), (31) 

y { ^r 2 + Az 2 J y 

where 

11 



_ yul x e y e m' h + x h y h m' e , 

A — — ; ; : — , /ill 



m' e m' h 



(32) 



y(m' e + m' h ) ' 11 m' e + 

Clearly, e(A) is (in units of Ry') the ground-state energy of an anisotropic Coulomb system. 
Before we comment on this, we turn to the last term in inequality (p9|), namely e' im . We 
transform the corresponding operator h' im in two steps. Firstly, we replace Zh by —Zh, 
secondly, we introduce center-of-mass and relative coordinates instead of z e and Zh as was 
done above. We find 



inf spec 



x e (l - y e ) x h (l- y h ) 



m' 



my, 



dz 2 



2q 



1 1 

R— R-i 



, (33) 



where R + and i?_ are derived from r + and r_ by replacing z e + Zh by the relative coordinate 
z. Splitting the right-hand side of expression (|33"D again into two parts, containing the terms 
and 1/R+ separately, we may shift the variable z by +1' and by —V, respectively. We 
obtain pure Coulomb potentials in both cases. Assuming their solutions to be known, we 
are finally lead to the inequality 



e' im > inf spec <^ -(1 - y)V| - 



XeQ- ~ Ve) X h (l - y h ) 



m' 



m. 



d 2 Aq' 

dz 2 r 



(34) 



Performing an appropriate scaling of variables, we can reduce the right-hand side once 
more to the anisotropic Coulomb problem. In comparison with Eq. ( |3i~D , the parameters are 
changed: 



e- > 

ira 



4q 2 



-inf spec < —V 2 — 



1-2/ 

the anisotropy parameter now being 

li' x x e (l - y e )m h + x h (l - y h )m. 



rj + A'z 2 



(35) 



.4' 



(36) 



/ijl (1 - y)(m e + m h ) 

We summarize our results for e' c , e' im , and E /Ry'. To evaluate the relations ( p9|) , (|3T|) 
and (^), we need an exact expression for the ground-state eigenvalue e(A) of the anisotropic 
Coulomb problem or, at least, a corresponding lower bound. Assuming e(A) to be known, 
the above considerations prove the following inequality: 
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E . 

> max 

Ry' Xe,x h ,y e ,y h ,y 



e e + e' h + -e(A) + j^-MA') 



(37) 



Finally, we arrive at a result for the binding energy, which complements the relation (|2 



Ry' 



< 



max 

Xe,Xh,Ve,yh,y 



e , e ) + (e' h - e 0lh ) + -e{A) + -^—e(A') 



(38) 



To the best of our knowledge, an exact analytical equation for e(A) is not available up to 
now. However, involved variational calculations were performed by several authors; we refer 
to the paper by Gerlach and PollmannEl, the references therein and to our brief summary 
in the Appendix. The optimal upper bound on e(A), which is presented in Ref. [15], will be 
denoted as ecp(A). Numerical studies indicate that ecp(A) deviates from e(A) only on a 
one-percent scale, although a rigorous estimate is admittedly missing. We mention that the 
quoted paper also contains additional lower bounds for e(A) which might be used to evaluate 
the preceeding inequality. In this paper we shall add another lower bound (again, see the 
Appendix). Unfortunately, the overall quality of all these lower bounds is not sufficient. 
Therefore, we will usually insert the numerical approximation ecp(A) instead of e(A). 



IV. RESULTS AND DISCUSSION 
A. Absence of mass- and dielectric mismatch 

We used this simplified model to test the efficiency of the above bounds and to illustrate 
the dimensional aspect. The system is interesting on it's own and was already discussed 
in the literature (see Refs. 1,0). Here, it will prove useful to understand the well-width 
dependence of the binding energy on a quantitative basis. 

One can easily specify the general results from above for the discussion of the present 
case. Equalizing all primed and unprimed material parameters, we have q = and h = h'. 
We simplify even further by assuming m e = m h , fi± = m = m e /2, and U' e = U' h = U. 
Turning to inequality ( p8|) and recalling its derivation, we have to choose y e = = y = 1 
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(no image potential). Consequently, the anisotropy parameter (see Eq. (j32|)) is fixed as 
A = (x e + Xh)/2. Summarizing at this stage, we derive from relation (|38|): 
E b 

— < - max [(e e - e , e ) + (e h - e ,h) + e(A)} , (39) 

Jriy %e,Xh 

where ej (see Eq. ( p0[ ) is now the solution of the equation 

2(1 - x { 



p arcsin \Jl - d/ll, (40) 

and e ij0 can be derived from e^, if we let Xj = 0. The eigenvalue e(A) was defined in Eq. ([JT]). 
We shall now treat the cases of an infinite and a finite well. They will be shown to have a 
qualitatively different behavior in the narrow-well regime, as indicated in the Introduction. 
For the limiting case of an infinite well we can simplify the inequality (|3~9"D as follows: 

E b 



< — max 
Ry~ A 



7T 2 



(41) 



P 

To proceed, we employ another inequality, namely e(A) > —4/(1 + 3A), which is derived in 
the Appendix. In this case, we are lead to an analytical result for Ej,: 



Ry ~ 



47^ 



4> if I < ^ « 0.91, 
6 



r 2 



31 3/ 2 



"" if 0.91 </< 3.63, (42) 



7T 2 .. . . 2tt\/3 



1 + — , if I > w 3.63. 

The right-hand side of Eq. fl42|) is shown as dashed-dotted curve in Fig. [| One realizes 
that the required two- and three-dimensional limits 4 Ry and 1 Ry are reproduced, if I tends 
to and oo, respectively. Unfortunately, the overall quality of this bound is rather poor. 
Inserting the result ecp(A) for e(A) as indicated in the preceding section, the upper bound 
on Eb is drastically lowered as shown in the figure (dashed line for V = oo). In addition, 
Fig. |3] contains a lower bound for the binding energy (solid curve for V = oo), which is based 
on inequality (|25|). 

Before we discuss the results in greater detail, we consider the finite-well case. Then, the 
inequality (^) has to be evaluated without further simplifications. Apart from the limiting 

14 



cases I — > and I — > oo, an analytical discussion is not possible. The small-width case, 
however, is interesting on its own. For I — > 0, the solution for e$ is 

e; = C/[l - f// 2 /(2 - 2^) + 0(/ 3 )]. (43) 

Inserting this expression in inequality and utilizing e(A) > — 4/(1 + 3 A) again, one 
finds the result 

1 < ^ < 1 + VsUl + OU 2 ), if Z->0, (44) 
Ry 

which was mentioned in the abstract as well as the Introduction. For finite values of I, we 
proceed as in the numerical treatment of the infinite well. Replacing e(A) by ecp{A), we 
find the upper bounds, which are shown in Fig. | (dashed curves). Again, the figure includes 
the corresponding lower bounds (solid curves). 

Clearly, the most significant attribute of these curves is the appearance of a maximum of 
Ej, as function of the well width L, if the height V of the well is finite. There exists a unique 
relation between V and the position L max of the maximum; for larger values of V L max 
becomes smaller. All curves exhibit the expected asymptotic behavior. For finite (infinite) 
V, the large-L limit of the binding energy is lRy, the small-L limit 1 Ry (4Ry). Last but 
not least, the shape of the upper and lower bound is very similar, the relative deviation not 
exceeding 0.2. We are sure that the deviation as such is mostly due to the inaccuracy of 
the upper bound on the binding energy. We recall that the lower bound for Ej, is based on 
a variational upper bound on the ground-state energy, whereas the upper bound needs an 
accurate lower bound for E as input; as usual, this part of the task is the more difficult 
one. 

B. Finite rectangular quantum well in a heterostructure 

In this part we apply our theory to a single-well structure, which exhibits a mismatch of 
both masses and dielectric parameters. As far as specific material data are concerned, we 
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use those of Ga\- x Al x AsjGaAsjGa\- x Al x As. We stress, however, that a direct compari- 
son of our results with experimental data is limited by the fact that the present theory is 
clearly incomplete. An obvious shortcoming is that, for example, the effects of valence-band 
degeneracy, spin-orbit coupling, and exchange interaction are not included. Our intention 
was to analyze the implications of the well structure, in particular the inability to sepa- 
rate the center-of-mass and relative coordinates, as accurate as possible in order to have a 
well-defined basis for further improvements. 

For the material parameters, we refer to the tables of Landolt-Bornstein@ and the work 
of WinkleiH (see also references therein). For GaAs, we used e = 12.53, m e = 0.067mo, 
furthermore = 0.090m , = 0.051m , Ry = 4.418 meV, a B — 130.21 A for the 
light hole, and m h = 0.377m , = 0.042m , Ry = 3.638 meV, a B = 158.12 A for the 
heavy hole. For AlAs, the corresponding parameters are e = 10.06 , m e = 0.150m , = 
0.208 (0.478)m , fjt± = 0.106 (0.093)m , Ry = 14.2(12.5) meV, a B = 50.2(57.2) A for the 
light (heavy) hole. The material parameters for Ga\_ x Al x As are normally described by 
linear interpolation formulae, for instance, e' = 12.53(1 — x) + 10.06a;, with one notable 
exception; the difference AE g of the gap-energies for x = and x > is fitted asi 

AE g = (1.087a; + 0.438x 2 ) eV. (45) 

The precise partition of AE g on the electronic part V e and the hole part Vh has been 
critically discussed; to the best of our knowledge, no previous general solution exists. We 
used the empirical relation V e = 0.65 AE g and Vh = 0.35 AE g (again, see Ref. |9|), but also 
indicate the consequences of modification to this assumption. 

We discuss now our results. To begin with, we present lower bounds for the binding 
energy according to relation (p5|). In Fig. [| we depict the influence of the parameter mis- 
matches on the heavy-hole (HH) exciton of Al x Ga\_ x As j GaAs j Al x Ga\_ x As for x = 0.3. 
The thin solid curve describes the hypothetical case of equal masses and dielectric constants 
in well and barrier; it may be viewed as a reference line. Switching on the discontinuity of 
the band masses only (circles), the bound for the binding energy is lowered. We are confi- 
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dent that this is an artifact of our variational treatment. On the one hand, the energy of 
the continuum edge is exact, the effects of the mass mismatch being fully incorporated; on 
the other hand, the variational bound on the ground-state energy is, of course, an approxi- 
mation. The difference between the two underestimates the influence of the mass mismatch. 
Considering only a dielectric mismatch, we find the curve depicted by triangles. One realizes 
that the presence of image charges shifts the binding energy to higher values, and this shift 
remains present for a wider range of the well width. It is only for L ^> as that the energy 
shift disappears. Finally, the thick solid line summarizes all effects. In comparison to the 
reference curve, the peak height and peak position is only slightly changed; in our case, we 
find L max ~ 60 A and a peak height of 9.5 meV. 

Fig. |5] presents lower bounds for the binding energies of heavy- and light-hole excitons 
in Al x Gai_ x As/GaAs/Al x Gai_ x As, assuming different ^/-concentrations x = 0.15, x = 
0.3, x = 0.4. The dashed (solid) curves correspond to the light-hole (heavy-hole) exciton. 
The lowest (highest) curve of every set belongs to the smallest (largest) value of x. The 
observed sequence can easily be explained. A higher ^/-concentration causes a higher band- 
gap difference and therefore higher potential wells. In addition, the peak positions are 
found to decrease with increasing ^/-concentration. Again, the reason is clear; the higher 
the value of x, the smaller the influence of the barrier. In fact, we expect to recover the 
trend illustrated in Fig. |3|, which is indeed the case. Above all, we find a peak structure for 
x > 0.25 and L > 30 A, in contrast to the assertion of Andreani and Pasquarelloll. 

We mentioned above that the actual values of V e and Vh are somewhat controversial. 
Therefore, we found it interesting to change the fraction Vh/V e . In Fig. |6| we compare lower 
bounds for the binding energy, Vh/V e being chosen as 35/65 (circles) and 15/85 (triangles), 
respectively; the latter value was used by Greene, Bajaj and PhelpsE Interestingly enough, 
the peak position is not changed too much. 

We shall now comment on the accuracy of our results. To do so, we have chosen the HH- 
case of Al 3 Ga 7 As/GaAs/Al 3 Gao, 7 As as a representative example. In Fig. [7| we contrast 
lower and upper bounds for the binding energy according to Eqs. ( P5]) and ([38]). The shape 
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of the curves and the peak positions are nearly the same for both bounds, the peak heights, 
however, differ disappointingly by a factor 1.5. Analyzing the numerical data in detail, we 
realized two points: i) the adapted trial function is well localized within the well for all 
values of L > L max and even slightly below; ii) the upper bound for the binding energy 
(or, equivalently, the lower bound for the ground-state energy) grossly overestimates the 
mismatch. If we take (i) for granted, also for the exact wave function, we can simplify 
the peak problem from the very beginning. If we calculate an expectation value of h (see 
Eq. fll4]) ) for such a wave function, we can firstly omit all 0-factors, and secondly replace all 
material parameters by well parameters. Then, h may be replaced by hf (see Eq. (f27D) with 
the important modification that all parameters in h! have to be understood as unprimed 
ones, namely those of the well. Calculating lower bounds in exactly the same manner as 
before, we find the bounds which are denoted by approximate ones (triangles). Of course, 
this curve has to be cut in the region of small well widths when the requirement of prevailing 
localization inside the well is not fulfilled any longer. One observes that the allowed channel 
for the exact binding energy is now considerably smaller. 

Finally, Figs. [5], [] compare experimental and theoretical results, again for 
Ga , 7 Al 3 As/ GaAs/ Ga ,jAl 3 As. The observed peak structure is reproduced by the present 
theory, but the experimental binding energies are up to 1 meV larger. There are indications 
that this discrepancy may be caused by our simplification of the real band structure. An 
enhancement of the exciton binding energies due to valence band degeneracy and conduction 



band nonparabolicity has been found by several authors (see, e.g., Refs. |3|j20D. 



V. CONCLUSIONS 

The intention of this work was to calculate the binding energy of an exciton in a rectan- 
gular quantum well with a higher reliability as was previously to be found in the literature. 
In doing so, we provided upper and lower bounds for the binding energy, which have a 
similar shape and constitute an allowed channel for the exact binding energy which is satis- 
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factorily small. We can conclude that the peaks of the binding energy as function of the well 
width are well established within our model, and for a wide parameter range. In the case 
of Gai- x Al x As/GaAs/Ga\- x Al x As, our results revise those of Ref. |8| thereby extending the 
regime of "admissable" well widths to less than 30 A. From a methodological point of view, 
our derivation of upper bounds on the binding energy (or, equivalently, lower bounds to the 
ground-state energy) can be transferred to related problems. 
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APPENDIX: UPPER AND LOWER BOUNDS FOR THE ANISOTROPIC 

COULOMB PROBLEM 

To evaluate the expression (|38D for the lower bound on the binding energy of Hamiltonian 
]), we had to provide an expression for the ground-state energy e(A) of the Hamiltonian 

2 



h A := -V 



2 



r 2 + Az 2 
2 



-V 2 - . (Al) 
lAr 2 + (1 - A)r\ 



We shall now briefly comment on this problem. Clearly, e(A) interpolates between the two- 
and three-dimensional limits of the hydrogen system. If A is changed from 1 to 0, the 
ground-state energy varies from -1 to -4 (in Rydberg units). Therefore, it is tempting 



to choose a variational wave function for the ground state as follows (see Ref. [151) : 

$ = Ce- A(r+ar±) , (A2) 
19 



where a, A are variational parameters. This wave function is asymptotically exact in the 
three- (A — > 1,A — > 1, a — > 0) and two-dimensional (A — > 0, Aa — > 2, A — > 0) limits. 
Calculating the expectation value of Hamiltonian ( |A1| ), one finds e(v4) < ecp(^4), where 

^ 2 (A,a) 



and 



e G p(A) = min(a 



tt/2 

V(x,a) = f d6- 
o 1 



V 2 (l,a) -a 2 \/ 2 (0,a)/4' 



sin# 



(A3) 



(A4) 



'1 + a sin 9) 2 V sin 2 # + x cos 2 
For a derivation of these equations and a discussion of the quality of the bound ecp(A), 
we again refer to Ref. [TBI. In Sec. IV we made repeated use of formula ( |A3| ). 

Interestingly enough, the same trial wave function provides us with a lower bound for 
e(A). To demonstrate this, we start from the identity 



V 2 ^ + V^ 

V tr 



-A 2 (l + a) 2 V, 



2A a\ , 9 / r I 
— + — + 2aA 2 1 - — 
r rj_ V r 



(A5) 



Thus, ?/> is an eigenfunction of the Hamiltonian h tr = —V 2 + V tr , corresponding to the 
eigenvalue e tr = — A 2 (l + a) 2 . In fact, ip is the ground-state eigenfunction of h tr . To prove 
this, one should realize that i) the ground-state of ht r is nondegenerate for all a, and ii) ej r 
coincides with the ground-state energy for a = 0; besides, the chosen if) has no zeros. We 
remark that h tr has the same two- and three-dimensional limits as Ha- 

Now let tpA be the exact ground-state wave function of the Hamiltonian Ha- Then we 
obtain 



e(A) = (V>aMVu) 

= (lpA\h t r ~ 



^Ar 2 + (1 - A)r\ 



- V tr \lpA) 



. , ,2A aA _ , 2 /., r_L 
> e tr + — + — + 2aA 2 1 - — 
r rj_ V r 



> e tr + (iPa\~ 
r 



2A + ^ 
P 



A + (1 — A)p 2 



Ar 2 + (1 - A)ri 
|V>a>, 



(A6) 
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where we defined p := r±/r and made use of p < 1. The right-hand side of the latter 
inequality ( |A6[ ) may be considered as a function of the variational parameters A and aX 
(A being the experimental "input"). We evaluate this function as follows. Firstly, we 
calculate the minimum of the expression in square brackets as a function of p; inserting the 
corresponding solution into inequality ( |A6[ ), the square bracket is a pure c-number and can 
be extracted from the expectation value. Secondly, we choose A such that the extracted 
square bracket vanishes and, finally, aX such that the lower bound assumes a maximum. 
This leads us to 

e(A) > (A7) 

Obviously, this lower bound gives the correct values e(l) = —1 and e(0) = —4. 
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FIG. 1. Geometry of the quantum well 
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FIG. 2. Upper (solid lines) and lower (dashed lines) bounds for the anisotropic Coulomb prob- 
lem. The solution which is considered as numerically exact is drawn as a thick line. 
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FIG. 3. Binding energy as function of the well thickness for different potential heights, showing 
comparison of variational calculation (solid lines) and lower bound method (dashed lines). The 
dashed-dotted line shows the result for the simple analytical lower bound from the Appendix. 
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FIG. 4. Binding energy of a HH-exciton in Alo.3Gao.7As/GaAs/Alo.3Gao.7As; the effect of 
dielectric and mass mismatch is shown separately by the triangle and circle symbols, respectively. 
The thick solid line illustrates the result of both contributions. 
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FIG. 5. Binding energies of the HH- and LH-exciton for different ^-concentrations 
x = 0.15, x = 0.3, x = 0.4. The variational approach including parameter mismatch is applied. 
The lowest (highest) curves belong to the smallest (largest) value of x. 
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FIG. 6. Binding energies of the HH- and LH-exciton for different valence band offsets in 
AIqaGclq.qAs / GaAs / AIqaGclq.qAs. The valence band offsets are chosen either from Greene, Bajaj 
and PhelpJ (triangles) or from Andreani and Pasquarelloi (circles). 
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FIG. 7. Rigorous upper bounds (circles) and lower bounds (lines) in comparison with approxi- 
mate upper bounds for the exciton binding energy in Alo.3Gao.7As/GaAs/Alo.3Gao.7As. 
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FIG. 8. Comparison of the present calculations (solid and dashed-dotted lines) with exper- 
imental data (stars and triangles) as well as with theoretical data (square symbols) for the 
heavy-hole-exciton in Alo^GaojAs/GaAs/AlQ^Gao^As. 
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FIG. 9. Comparison of the present calculations (solid and dashed-dotted lines) with experi- 
mental data (stars) as well as with theoretical data (square symbols) for the light-hole-exciton in 
Alo.3Gao.7As/GaAs/Alo.3Gao.rAs. 
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